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10 WE: 涡 致 混合 扩散 是 物理 海洋 研究 中 的 热点 和 难点 问题 。 本 文 基于 “有 效 扩散 ”理论 ， 研 究 示 踪 物 等 值 
11 线 在 海 表 地 转 消 流 的 多 尺度 搅拌 作用 下 ， 发 生 拉 伸 、 扭 曲 、 变 形 、 折 肢 等 改变 其 几何 拓扑 结构 的 现象 ， 并 
12 ”探讨 了 等 值 线 分 形 长 度 的 变化 与 混合 效率 的 关系 。 研 究 结 果 表 明 ， 在 地 转 流 场 的 搅拌 下 ， 示 踪 物 的 等 值 线 
13 ”会 被 迅速 拉 长 ， 并 产生 丰富 的 精细 结构 。 这 种 分 形式 的 增长 可 达 原 长 度 的 10-20 倍 ， 是 混合 效率 提高 的 主 
14 ”要 原因 ; 而 涡 丝 和 锋面 伴随 的 梯度 增强 虽然 也 有 贡献 ， 但 为 次 要 因素 。 另 一 方面 ， 在 示 踪 物 模拟 过 程 中 
15 小 尺度 扩散 会 通过 不 可 逆 混 合 对 示 踪 物 进行 均匀 化 ， 从 而 抹 平 等 值 线 的 精细 结构 ， 抑 制 等 值 线 的 增长 ， 限 
判 混合 效率 的 提高 。 基 于 '“ 数 盒子 "算法 计算 了 等 值 线 的 分 形 维度 ， 其 数值 在 1.4 到 1.6 之 间 ， 介 于 一 维和 二 
17 ， 维 之 间 。 但 由 于 地 转 消 流 数据 分 辨 率 的 限制 ， 无 法 考虑 更 小 尺度 《次 中 尺度 过 程 ) 的 搅拌 作用 ， 可 能 低估 
18 ”了 等 值 线 的 分 形 长 度 和 混合 效率 。 本 研究 将 海洋 混合 与 等 值 线 几 何 特征 联系 了 起 来 ， 初 步 得 到 了 分 形 长 度 
19 ”和 混合 效率 两 者 的 经 验 关系 式 ， 未 来 可 以 利用 图 像 识 别 等 成 熟 遥 感 技术 将 海洋 示 踪 物 等 值 线 的 几何 特征 直 
20 ” 接 转 换 为 混合 效率 ， 为 诊断 分 析 海 洋 混合 及 其 参数 化 提供 了 一 种 新 的 思路 。 
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Abstract: Quantifying eddy mixing in the ocean is a hot and tough problem in the area of physical oceanography. 
Based on the theory of effective diffusivity, the present study investigated the stirring effects of geostrophic 
turbulence that lead to stretching, distorting, deforming, and folding of tracer contours. These changes are then 
related to the efficiency of turbulent mixing. Results show that under the stirring effect of geostrophic turbulence, 
the length of tracer contour can be quickly elongated and fine-scale tracer filaments and fronts are also generated. 
This fractal elongation of tracer contour, about 10-20 times longer than the original length, is the dominant 
contributor to the mixing efficiency, whereas the gradient enhancement associated with filament and front 
generations only plays a secondary role. On the other hand, fine-scale features are smoothed out by small-scale 
diffusivity which eventually suppresses the increase of contour length and the generation of tracer filaments. This 
imposes an upper bound of the mixing efficiency when the stirring and smoothing effects are in a dynamical balance. 
Through a ‘box-counting’ method, the fractal dimension of tracer contour is also found between 1.4~1.6, indicating 
a geometric dimension lies somewhere between 1D and 2D. Due to the limitation of data resolution, contour length 
and thus mixing efficiency may be underestimated. Finally, the present study made an empirical relation between 
the fractal dimension and mixing efficiency, providing an opportunity for estimating mixing efficiency through a 
well-developed pattern recognition technique in remote sensing, and a new way of diagnosing ocean mixing and its 


parameterization. 
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涡 致 混合 扩散 是 海洋 
动量 和 热量 , 通 


到 当前 的 气候 状态 。 研 究 指 出 ,混合 


物理 过 程 (Sallée et al., 2008), FEY 


此 外 ， 混 


合 扩散 系数 是 数值 模型 的 重要 


合 扩散 不 仅 是 引起 海洋 内 部 水 团 变 性 、 模 态 水 形成 的 重要 


过 程 . 从 整体 上 看 ,海洋 从 上 表面 获得 不 均匀 的 物质 、 


过 海洋 中 不 同 尺度 涡 旋 的 搅拌 运动 ， 使 其 在 各 个 方向 上 进行 高 效 混合 ， 以 达 


由 
pa 


洋 经 向 翻转 环流 的 关键 因素 (Marshall et al., 2017). 


寻找 最 优 或 合适 的 数值 ， 即 混合 扩散 参数 化 ， 


是 提高 模拟 结果 的 重要 步骤 (Monti and Leuzzi, 2010; Haza et al., 2012; Fox-Kemper et al., 2014; 


邢 元 明 等 , 2013 )。 


al., 2022), Wl 


此 也 是 物理 海洋 学 


合 扩散 具有 复杂 的 时 空 


分 布 特征 (Kamenkovich et al., 2021; Qian et 


海洋 中 的 运动 涵盖 了 极其 广泛 的 时 空 尺 度 。 地 球 的 大 小 限制 了 海洋 运动 的 最 大 空间 尺度 


( 约 10m)， 即 行星 尺度 运动 ， 而 分 子 间 的 粘性 则 决定 了 海洋 运动 的 最 小 尺度 〈 约 103m)， 


两 者 之 间 跨 度 达 


潮流 的 结构 特色 


10 个 数量 级 。 


满 流 作为 一 种 多 尺度 运动 ， 能 量 在 不 同 尺 度 上 的 分 布 决 定 了 


E， 特 别 是 在 中 尺度 (A 50~100km) 以 下 ， 运 动 具有 一 定 的 自 相 似 性 (Fox- 


Kemper and Menemenlis, 2008 )。 正 如 数值 天 气 预 报 的 先驱 之 一 Lewis Richardson F 1922 年 


写 的 诗 所 


性 耗 散 ”。 


述 的 那样 : 


正 是 这 种 自 相 似 性 


包含 小 涡 ， 并 喂 予 速度 ， 


小 涡 包 含 更 小 的 涡 ， 如 此 继续 直到 寿 


流 运动 与 分 形 几何 联系 起 来 。 


上 世纪 60 年 代 ，Mandelbrot (1967) 提出 了 分 形 几 何 学 来 描述 和 研究 表面 粗糙 、 具 有 自 


相似 性 质 的 几何 图 


形 。 他 提出 了 经 典 问题 “ 英 


的 海岸 线 有 多 长 ”并 指出 海岸 线 的 长 度 依 


赖 于 用 于 测量 的 太子， 尺子 越 短 ,测量 结果 越 长 。 这 


岸 线 放大 后 ,仍然 会 出 现 和 原来 相似 的 细节 ,。 因 


是 因为 海岸 线 也 有 自 相似 性 : 粗糙 的 海 


此 , 较 短 的 尺子 可 以 测量 出 更 多 细节 的 长 度 。 
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1 示 踪 物 等 值 线 长 度 与 混合 效率 关系 示意 图 。 

(a) 等 值 线 为 球面 上 最 短 曲线 时 的 分 布 (为 同心 圆 或 纬 圈 线 )， 此 时 混合 效率 最 低 ;(b) 气候 态 〈 欧 拉 时 间 平 均 ) 等 值 
线 分 布 ; 〈c) 瞬时 等 值 线 分 布 图 ， 此 时 其 混合 效率 相 比 (a) 和 (b) 具有 数量 级 的 提高 。 三 个 状态 中 等 值 线 包 围 的 面积 
相等 。 审 图 号 为 GS(2016)1665 


qi 


Fig.1 Schematic illustration of the relation between contour length of a tracer and mixing efficiency. 
(a) minimum possible contour length on a sphere (latitude circles) which is at the lowest mixing efficiency; (b) smooth 
climatological contour through a Eulerian time mean; (c) instantaneous wavy contour distribution in which mixing efficiency is 


greatly enhanced. Note that the areas between any two contours are the same for each panel. 


RE Mandelbrot 早 在 1975 年 就 注意 到 了 灌流 混合 与 分 形 几何 的 联系 (参见 黄 真理 , 2000; 


` 


Li et al., 2001; 刘 式 达 等 , 2014; Paladin and Vulpiani, 1986; 沈 学 会 等 , 2005 )， 在 海洋 领域 也 


有 将 分 形 几 何 应 用 到 漂流 浮标 轨迹 (Osbome et al., 1989; Sanderson and Booth, 1991)、 海 面 


Se, 


波浪 〈 田 纪 伟 等 , 1996; HAE, 2015) 等 海洋 环境 要 素 的 分 析 中 《〈 付 垦 华 ,2000)， 但 相关 研究 


仍然 有 限 。 此 外 , 还 没有 理论 将 混合 扩散 系数 、 混 合 效率 ( 涡 扩 散 系 数 与 分 子 扩散 系数 之 比 ) 


等 概念 与 等 值 线 的 几何 特征 联系 起 来 ， 仅 有 的 一 些 研 究 常 常 聚 焦 于 均匀 或 各 向 同性 的 消 流 


(Iyeretal., 2020)， 与 真实 海洋 的 复杂 情况 有 很 大 不 同 。Nakamura (1996) 提出 了 “有 效 扩 


散 ” 理 论 , 将 大 气 海洋 中 的 侧 向 混合 扩散 系数 与 等 值 线 的 拓扑 形状 联系 了 起 来 ,根据 该 理论 ， 
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22 


SP) iMi WE foe) WISH A EAI CAI la)， 不 同性 质 流体 的 


接触 面 ( 即 等 值 线 长 度 ) 最 小 ,扩散 效率 最 低 ; 在 粗 分 辨 率 海洋 气候 模式 中 ， 等 值 线 常 常 表 


现 为 一 组 光滑 的 曲线 (图 1b); 而 从 实际 的 高 分 辩 率 图 像 来 看 ， 等 值 线 并 不 光滑 ， 瞬 时 图 像 


十 分 复杂 。 当 等 值 线 被 流 场 搅拌 时 , 不 同 尺度 的 涡 旋 会 导致 等 值 线 分 形式 地 拉 长 ， 产生 丰富 


的 精细 结构 ， 如 锋面 和 涡 丝 、 湛 流 界 面 的 卷 衫 、 破 碎 且 细小 的 “孤岛 ”或 “气泡 ”( 图 1c)。 


通过 这 种 方式 极 大 地 扩大 了 不 同性 质 流体 的 接触 面积 ， 从 而 有 效 提升 了 混合 效率 。 


在 数值 模式 中 ,衡量 等 值 线 长 度 的 最 小 尺子 是 模式 网 格 的 边 长 ， 即 分 辨 率 。 对 于 低 分 辨 


率 的 气候 模式 ， 尺 子 较 长 ， 因 而 测量 到 的 等 值 线 具 有 简单 、 平 滑 的 几何 特征 ， 随 着 分 辨 率 的 


逐渐 提升 ， 尺 子 越 来 越 短 ， 从 而 能 扑 欣 更 多 丰富 的 等 值 线 拓扑 细节 。 对 于 流体 的 模拟 ， 存 在 


尺子 长 度 的 下 限 ， 也 就 是 分 子 粘 性 或 扩散 尺度 (Batcherlor 尺度 )。 小 于 这 个 尺度 的 运动 引发 


的 等 值 线 拓 扑 细节 会 被 粘性 或 扩散 迅速 抹 平 ， 因 而 在 这 个 尺度 下 ， 难 以 出 现 更 丰富 的 细节 。 


正 是 因为 低 分 辨 率 模式 中 的 混合 效率 比 高 分 辨 率 模式 要 低 , 常常 需要 设置 更 大 的 扩散 系数 来 


进行 耗 散 (Huang, 2014). 


本 研究 尝试 将 数学 领域 的 分 形 几 何 理论 (Mandelbrot, 1967) 和 地 球 流体 领域 的 “有 效 扩 


BY” 到 


EY (Nakamura, 1996) 联系 起 来 ， 利 用 卫星 观测 的 海 表 地 转 流 场 驱动 被 动 示 踪 物 ， 通 


过 示 踪 物 的 等 值 线 信息 , 考察 全 球 海 表 涡 致 混合 的 效率 和 等 值 线 几何 特征 之 间 的 关系 ,以 期 


给 出 两 者 之 间 的 可 能 联系 本 文 的 第 一 部 分 介绍 数据 和 分 析 方 法 ,第 二 部 分 是 诊断 分 析 结 果 ， 


最 后 是 结论 和 讨论 。 


1 数 


据 与 分 析 方 法 


1.1 有 效 扩散 理论 


4 


FE 意 一 个 被 动 示 踪 物 4 在 流 场 u ASK Pie ALERT AE ACE CA, o ) 平流 -扩散 方程 : 
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<4(t,4,9) = Vg+t rv’g. (1) 


其 中 ，u = (u,v) 是 水 平流 场 ，4 EAR, p 是 纬度 ，K, 是 次 网 格 扩 散 系 数 ，V 是 球面 二 维 


梯度 算 子 。Nakamura (1996) 提出 用 示 踪 物 的 一 组 等 值 线 g 包围 的 面积 


ala) = ff ,dS, (2) 


作为 新 的 坐标 来 描述 示 踪 物 的 运动 。 公 式 (2) F dS =a’ cosedAdg 是 球面 面积 元 ， 


q (4,0) 为 某 时 刻 示 踪 物 的 二 维 球面 分 布 场 ，a 为 地 球 半径 。 定 义 等 值 线 积分 算 子 


(= 外。 ()4S， 风 有 4(9) = 4(1) 。 由 于 面积 坐标 不 符合 通常 的 习惯 ，Haynes 和 


Shuckburgh (2000) 进一步 提出 了 相当 纬度 p。 (4(9)) = Pu (9g) 来 蔡 代 面 积 坐 标 .4。 相 当 续 


度 的 物理 意义 是 , 该 纬度 线 以 南 包 围 的 海水 面积 与 对 应 的 等 值 线 包 围 的 海水 面积 相等 。 这 样 
相当 纬度 与 等 值 线 就 满足 一 一 映射 的 关系 , 可 以 用 来 做 空间 坐标 。 他 们 进一步 给 出 了 平流 扩 
散 方程 (1) 在 相当 纬度 坐标 下 的 形式 : 


l 
o K œ coso 9 (3) 
eff eq 
~ cos Pog apu a0 Pa 


方程 (3) 表明 ， 变 换 空间 坐标 到 相当 坐标 系 Aa Pa) 后 ， Oy, 表示 垂直 于 等 值 线 的 广 


向 ， 而 A, 表示 沿 等 值 线 方向 。 这 种 变换 等 价 于 将 示 踪 物 浓度 进行 南北 方向 的 《绝热 ) 排序 


(图 la)， 最 终 得 到 与 纬 圈 平 行 的 等 值 线 分 布 ， 并 且 单 调 递增 或 递减 《类似 于 垂 向 排序 得 到 


最 小 位 能 状态 )。 由 于 沿 等 值 线 方向 没有 变化 , O/ OA,, = 0 ,因此 二 维 空间 问题 简化 为 一 维 。 


此 外 ， 由 于 平流 作用 只 改变 等 值 线形 状 , 不 改变 等 值 线 包围 的 面积 ， 因 此 平流 项 不 再 显 式 出 


现在 方程 G) 中。 最 后 ， 扩 散 系 数 从 分 子 扩散 4 变 成 了 有 效 扩散 天 ss ， 其 定义 为 : 


=k eq (4) 


其 中 工 是 等 值 线 的 相当 长 度 ，ZLii, 是 最 小 可 能 长 度 ， 在 球面 上 即 纬 图 长 度 。 它 们 比值 的 平 


¢ 


Fi Leg | Ls ， 表 示 通 过 平流 搅拌 使 得 分 子 扩散 系数 被 放大 的 倍数 KK / Kp HU AHI BR 


[x 


Fl Fe 


M (Nakamura, 2001) 或 Nusselt 2 (Marshall etal.,2006)。 这 里 的 混合 效率 定义 


向 清流 混合 中 的 定义 ， 即 “动能 耗 散 向 势能 转化 的 比例 ” 当 等 值 线 分 布 与 图 la 一 致 时 ， 
6 


22 


Li, = 过， 因此 有 效 扩散 简化 为 分 子 扩散 ， 混 合 效率 最 低 (为 1 );， 当 等 值 线 被 灌流 充分 搅 


fia» L > L 


‘min ” 


此 时 有 kK, >K, ， 混 合 效率 可 以 使 有 效 扩散 系数 比分 子 扩散 系数 大 若 


干 个 数量 级 。 


相当 长 度 可 以 根据 下 式 (Nakamura, 1996; Haynes and Shuckburgh, 2000) 计算 得 到 : 


GAY Ii aay 8 f 
EE o 


其 中 (), =0A(-)/0A 是 g 临近 两 根 等 值 线 间 的 (宽度 加 权 ) 平均 算 子 。 最 小 可 能 长 度 志 


min 


则 由 排序 后 等 值 线 所 在 纬度 的 长 度 决 定 。 由 于 有 陆地 存在 , 纬 圈 长 度 由 所 有 网 格 中 的 海水 网 


格 宽度 累加 得 到 : 


Ln (q) = da COS Por, > (6) 


此 外 , 等 值 线 的 实际 长 度 志 =$., dl ,也 叫 等 值 线 周 长 , 可 以 通过 分 段 线性 累加 计算 得 到 ， 


具体 分 段 坐标 由 “行进 网 格 ”(marching squares) 算法 (Lorensen and Cline, 1987) 计算 获得 ， 


然后 将 两 点 间 的 坐标 转换 到 球面 距离 再 累加 起 来 即 可 算出 工 。 三 个 等 值 线 的 特征 长 度 满 足 不 


FAXR Ly Z L2 Lain (Haynes and Shuckburgh, 2000 )。 


实际 诊断 计算 步骤 简单 概括 如 下 : 首先 ， 对 于 任意 时 刻 示 踪 物 的 空间 分 布 ， 在 最 小 和 最 


大 值 之 间 选 取 一 组 离散 的 等 值 线 数值 g, ， 然 后 依次 计算 这 组 等 值 线 包围 的 水 点 面积 4(g, )、 


等 值 线 长 度 上 (gq,) 、 等 值 线 相当 长 度 L (q) 。 之 后 再 选取 一 组 纬度 等 值 线 g, ， 计 算 它们 包 


围 的 面积 4 (9,) ， 根 据 两 组 离散 值 4(9,) A A(Q,) 面积 相等 ， 可 以 通过 插值 得 到 离散 的 


Pa (4; ) 坐标 转换 关系 式 ， 最 后 利用 V。 (9; ) 关系 和 公式 〈6) WREE Lan (9 ) 。 本 文 主要 


考察 在 地 转 渍 流 搅 拌 下 ， 等 值 线 的 实际 长 度 工 与 相当 长 度 工 , 的 演变 情况 ， 进 而 探索 等 值 线 


分 形 拓扑 结构 与 混合 效率 到 / Li 的 关系 。 


1.2 地 转 流 场 数据 及 示 踪 物 离线 模拟 


本 文 利用 美国 床 省 理工 学 院 开发 的 大 气 -海洋 通用 环流 模式 MITgcm (MIT General 


7 


Circulation Model, Marshalletal., 1997) 对 示 踪 物 进 行 离线 (offline) 式 的 模拟 ， 从 而 得 到 示 


踪 物 等 值 线 随时 间 的 演变 数据 。 离 线 模拟 无 需 动力 内 核 ， 在 给 定 示 踪 物 初始 场 的 情况 下 ， 数 


值 求解 示 踪 物 控 制 方程 为 (1)。 而 驱动 示 踪 物 演变 的 海 表 流 场 来 自 于 AVISO 提供 的 网 格 化 


fal 


卫星 高 度 计 融合 产品 〈http:/www.aviso.altimetry.fr/duacs/)， 该 数据 集 包 括 空 间 分 辩 率 为 1/4 


度 、 时 间 分 辩 率 为 一 天 的 海 表 高 度 异 常 (SLA)、 绝 对 动力 地 形 (ADT) 以 及 根据 地 转 关 系 


计算 得 到 的 海 表 地 转 流 : 


u =(u sv )=4 -a n i (7) 
8 E g f Oy Ox 


其 中 ,是 地 转 流 ，7 是 卫星 观测 的 海 表 高 度 异常 ，g 是 重力 加 速度 ， 了 是 科 氏 参数 。 在 赤 


道 附近 fT 0， 该 产品 采用 Lagerloef (1999) 年 提出 的 方法 计算 赤道 附近 【〈 土 5° 范 围 


内 ) 的 地 转 流 。 为 了 减少 截断 误差 ， 本 文 首先 将 地 转 流 插 值 到 0.1” 分 辨 率 ， 与 模型 分 辩 率 


一 致 ， 采用 900s 积分 步 长 。 根据 Shuckburgh 等 (2009) 的 研究 , 海 表 流 场 符合 “紊乱 平流 ” 


特征 ， 因 此 插值 所 带 来 的 网 格 尺度 的 不 确定 性 不 会 严重 影响 本 文 的 计算 。 此 外 ， 由 于 地 转 流 


并 不 严格 满足 球面 无 辐 散 要 求 ， 本 文采 用 Marshall 等 (2006) 的 方法 ， 对 流 场 进 行 了 散 度 订 


正 ， 使 之 满足 无 辐 散 约束 。 因 此 ， 在 没有 源 和 汇 的 条 件 下 ， 示 踪 物 等 值 线 包 围 的 面积 4 是 严 


格 守恒 的 。 


模式 积分 需要 给 定 示 踪 物 初始 场 。Abernathey 和 Marshall (2013) 使 用 了 三 种 初始 场 ， 


分 别 是 纬 向 均匀 分 布 场 〈 图 1a)、 气 候 态 海 表 温度 分 布 场 、 气 候 态 流 函 数 分 布 场 。 三 种 初始 


场 中 ， 纬 向 均匀 分 布 场 更 适合 本 研究 的 问题 。 这 是 因为 它 直接 代表 了 最 低 混 合 效率 状态 ， 可 


以 直接 考察 在 地 转 消 流 搅拌 下 的 效率 逐渐 增加 的 过 程 。 其 次 , 纬 向 均匀 分 布 保证 了 示 踪 物 浓 


度 从 北 到 南 单调 增加 , 保证 了 相当 纬度 和 真实 纬度 之 间 有 对 应 的 统计 关系 , 这 样 在 相当 纬度 


空间 的 展示 也 基本 反映 了 真实 纬度 的 情况 。 基 于 上 述 考 虑 ， 本 文采 用 纬 向 均匀 分 布 作为 示 


踪 物 初始 场 〈 图 2a)。 
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1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 
图 2 ANREP BR 
初始 示 踪 物 浓 度 与 纬 圈 一 致 ， 且 被 标准 化 到 1-2 之 间 。 审 图 号 为 GS(2016)1665 


Fig.2 Evolution of a passive tracer driven by surface geostrophic flows. 


The tracer is initially aligned with latitudes and normalized between 1 and 2. 


图 2 AE SANE TR PRR TOL. EAEE i A EE Ps 大 约 20 天 后 (图 


2b)， 示 踪 物 就 产生 了 丰富 的 细 尺 度 结构 。 随 着 时 间 的 发 展 ， 在 海盆 尺 度 环流 的 作用 下 ， 原 


本 纬 向 分 布 的 等 值 线 (未 显示 标 出 ) 逐渐 开始 倾斜 : 北半球 的 等 值 线 汇聚 于 西边 界 的 湾流 区 


和 东边 界 的 热带 地 区 ,南半球 由 于 南极 绕 极 流 的 存在 ,等 值 线 只 是 在 低 纬 地 区 的 东边 界 汇聚 ， 


高 纬度 地 区 仍然 基本 维持 纬 向 分 布 。 此 外 ， 由 于 地 转 满 流 中 丰富 的 中 尺度 涡 旋 ,在 它们 的 捞 


拌 下 出 现 了 丰富 的 中 尺度 折 蕾 特征 , 使 得 等 值 线 不 再 光滑 平整 。 这 些 不 同 尺度 的 搅拌 作用 将 
会 极 大 地 提高 消 流 混合 的 效率 。 


1.3 分 形 维度 计算 方法 


等 值 线 的 分 形 维度 有 较 多 的 定义 , 不 同 定义 得 到 的 结果 也 存在 一 定 差异 。 本 文选 取 常 用 


1 “的 Minkowski 维度 来 表征 等 值 线 的 分 形 维度 ， 也 叫做 “ 数 盒 子 ”(box-counting) 维度 。 它 的 


2 ”计算 方法 采用 “ 数 盒子 ”法 ， 即 利用 一 定 分 辨 率 的 网 格 x 去 测量 等 值 线 ， 如 果 等 值 线 落 在 网 


3 ， 格 里 ,， 则 这 个 网 格 的 边 长 计 入 等 值 线 的 分 形 长 度 ， 否 则 不 计 入 ， 最 后 得 到 > 度量 下 的 盒子 数 


4 f(r) 以 及 等 值 线 长 度 L~x fr 。 根 据 Minkowski 分 形 维度 的 定义 〈 黄 真理 , 2000; Sreenivasan 


5 and Meneveau, 2006 ): 


è a __ log f(r) 


， (8) 
min log(r) 


7 ”需要 改变 网 格 分 辨 率 ”>， 即 用 一 组 不 同 长 度 的 太子 来 测量 等 值 线 长 度 ， 越 大 的 太子 则 会 忽略 


8 ” 越 多 的 细节 信息 。 通 过 拟 合 多 次 采样 得 到 的 等 值 线 长 度 和 分 辨 率 (或 者 尺子 ) 的 关系 ， 在 对 


9 ，” 数 坐标 图 上 得 到 曲线 斜率 ， 也 就 是 等 值 线 的 分 形 维度 Dj, 。 理 论 上 如 果 等 值 线 是 简单 光 清 


10 ”的 曲线 ， 则 分 形 维 度 为 1;， 如果 等 值 线 有 很 丰富 的 细节 被 越 来 越 小 的 盒子 捕捉 到 ， 则 分 形 维 


11 ”上 度 会 大 于 1， 但 不 会 超过 2， 即 介 于 一 维和 二 维 之 间 。 


12 实际 计算 中 ， 由 于 地 转 流 场 数据 插值 到 了 0.1”,“ 盒 子 ” 的 边 长 以 0.1° 为 基础 ， 逐 步 扩 


13 ”大 两 倍 ， 得 到 0.1*?、0.2?、0.4°*、0.8°、1.6*、3.2*、6.4*、12.8° 共 8 个 分 辨 紊 的 “盒子 ”。 这 


14 ” 样 可 以 在 对 数 图 上 得 到 8 个 等 间距 的 样本 ， 通 过 最 小 二 乘法 拟 合 得 到 Di, 。 通 过 多 次 计算 


15 ”表明 ， 在 最 大 和 最 小 分 辩 率 之 间 增 加 采样 分 状 率 并 不 会 显著 改变 拟 合 结果 。 


17 2 等 值 线 长 度 特征 


18 ”2.1 等 值 线 长 度 的 时 空 演变 特征 


Ww 
HÒ 


全 出 了 等 值 线 的 最 小 可 能 长 度 、 纬 圈 总 长 度 以 及 两 者 之 比 。 其 中 最 小 可 能 长 度 的 百 


19 øl 


20 ”分 比 (图 3b) 也 代表 了 海洋 网 格 与 总 纬 图 网 格 的 比值 。 海 洋 占 比 最 高 的 纬度 在 南极 绕 极 流 


21 ”区 域 ， 特别 是 Drake Passage 纬度 附近 (~50°S) 几乎 全 是 海洋 ， 在 赤道 附近 则 降低 到 75%. 
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从 20°N 开始 ， 比 例 逐 


减少 ， 且 其 纬 圈 占 比 也 同时 降低 ， 加 上 靠近 


Bit PRERIO 附近 。 特 别 是 在 南 、 北 


相当 纬度 (eq) 


-(b) — Linin/Liat 


0 20000 
KE /公里 
图 3 (a) AIK EL min (km, HR) 和 纬 圈 总 长 度 (km, R). Cb) Lmin AB 


70000 0 25 50 75 100 
百分比 /% 


EER 60? 以 外 ， 不 仅 海洋 网 格 迅速 


[极地 ， 示 踪 物 演变 受 陆 地 和 海 冰 影响 较 大 ,计算 


结果 有 一 定 不 确定 性 。 因 此 本 文 主要 关注 南北 纬 60° 以 内 的 计算 结果 。 


圈 总 长 度 的 百分比 。 


T 


Fig.3 (a) Minimum possible length Lmin (km, blue) and total length of latitudinal circle (km, orange), and (b) the 


ratio of Lmin to the total length of latitudinal circle. 


本 文 进 一 步 计 算 并 对 比 了 初始 时 刻 和 积分 一 年 后 等 值 


初始 时 刻 ， 等 值 线 为 纬 圈 线 。 以 图 4a 中 的 最 小 可 能 长 度 也 


线 的 长 度 。 从 图 4a 可 以 看 到 ， 在 


min 为 单位 ， 标 准 化 实际 长 度 和 相 


当 长 度 ， 那 么 两 者 都 是 1 (图 40)。 此 时 ， 混 合 效率 Mp = 已 / La 亦 为 最 低 值 1， 其 中 的 


微小 误差 并 不 显著 ,是 由 于 不 同 计算 方法 导致 的 。 在 积分 一 年 后 ， 等 值 线 在 各 个 斥 度 上 都 产 


生 了 丰富 的 细节 信息 


(图 4c)， 尤 其 是 在 赤道 东 太 平 洋 、 赤 道 印 度 洋 、 北 半球 西边 界 湾流 延 


展区 ， 由 于 不 稳定 波动 或 者 涡 旋 破 雄 的 搅拌 作用 ,等 值 线 断 和 裂 破碎 ,产生 了 局 部 复杂 的 拓扑 


结构 。 定 量 计算 表明 , 大 部 分 纬度 上 的 等 值 线 长 度 都 比 初始 时 刻 增长 了 10 倍 以 上 (图 4d)。 


尽管 整体 上 等 值 线 仍然 呈 东 西 走向 , 但 丰富 的 细节 充分 拓展 了 它们 的 长 度 。 其 中 南半球 副 热 


He (~30°S)、 赤 道 


(~3°N)、 北 半球 副热带 (~30°?N)、 北 


E 球 副 极 地 (~55°N) 有 四 个 局 部 


11 


峰值 (图 44)， 从 侧面 反映 了 南北 半球 涡 旋 活跃 的 区 域 以 及 赤道 波动 活路 区。 在 60°N 以 北 
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也 可 以 看 到 一 个 峰值 ， 这 可 能 与 格陵兰 岛 东 岸 的 急流 ( 东 格陵兰 流 ) AAR, MAAT 


陆地 和 海 冰 的 影响 也 存在 一 定 的 不 确定 性 。Bithell 和 Gray (1997) 对 大 气 对 流 层 顶 附近 的 


位 涡 等 值 线 也 做 了 类 似 的 分 析 。 他 们 的 计算 结果 表明 ， 对 流 层 顶 等 值 线 长 度 可 以 增加 约 10 


倍 以 上 ， 但 只 需要 一 周 时 间 。 


0° 60°E 120° 


60° 


4 示 踪 物 水 平分 布 及 等 值 线 长 度 


左 列 为 示 踪 物 的 水 平分 布 ， 黑 实 线 为 


E 意 选取 的 一 组 等 值 线 月 
12 


有 于 展示 ] 


相当 纬度 (hea) 


8 è è 


相当 纬度 (pea) 


SS k 


其 拓扑 细节 。 右 列 为 等 值 线 长 度 LC) 和 相当 


S è 8 


co 


(d) 5 


o è 8 o 8 8 8 


8 è 8 


o 


长 度 L。。( 检 色 )， 两 者 均 用 Li 做 单位 进行 无 量 纲 化 。 第 一 行为 初始 时 刻 的 情况 ， 后 三 行 分 别 为 用 小 尺度 扩散 系数 kc,, = 
5、20、50 ms! 积分 一 年 后 的 结果 。 左 边 图 中 审 图 号 为 GS(2016)1665 


Fig.4 Tracer horizontal distribution and its contour length L and equivalent length Leg. 


stg 


Left columns are tracer distributions with several contours highlighted. Right columns are contour lengths 
normalized by Lmin. The first row is the results at first day and the last three rows are results after 1-yr integration 


with small-scale diffusivity km 5, 20, 50 m's”. 


NO om FP WN FP 


8 对 比 等 值 线 的 实际 长 度 和 相当 长 度 , 可 以 看 到 在 初始 时 刻 (图 4b), 两 个 长 度 完全 一 致 。 


9 ”但 在 地 转 流 搅拌 下 (图 44)， 相 当 长 度 大 于 实际 长 度 ， 前 者 大 约 是 后 者 的 1.6 倍 ， 这 是 由 了 


10 ”大 量 涡 丝 和 锋面 相伴 随 的 梯度 增强 带 来 的 。2.2 小 节 将 进一步 分 析 两 者 之 间 的 关系 。 


14 上 面 看 到 的 等 值 线 的 拉 伸 和 增长 只 不 过 是 平流 -扩散 方程 (1 ) 中 平流 项 的 作用 ， 扩 散在 


12 ”等 值 线 的 演变 过 程 中 也 有 重要 作用 , 但 主要 在 网 格 尺 度 附近 。 因 此 ,在 示 踪 物产 生 丰 富 的 精 


13 ，” 细 结构 之 前 ， 难 以 看 出 扩散 的 作用 。 本 文采 用 不 同 的 小 尺度 扩散 系数 K, =5. 20. 50m*-s" 


14 ”进行 了 三 次 模拟 试验 。 图 4 展示 了 积分 一 年 后 示 踪 物 的 分 布 。 对 比 图 4c、4e、4g 可 以 看 到 ， 


15 ”KK 越 小 ， 则 等 值 线 越 多 精细 结构 ， 甚 至 在 部 分 区 域 缠绕 成 “意大利 面 ”的 形状 ， 当 x 逐渐 


m 


16 ” 增 大 时 ， 细 节 结 构 被 逐渐 抹 平 ， 到 了 为 50 ms 时 ， 等 值 线 的 精细 结构 所 剩 无 几 ， 但 仍 


17 ， 然 有 涡 旋 尺度 的 特征 。 从 等 值 线 的 长 度 也 不 难看 到 (图 4 AD, 扩散 系数 越 小 , 等 值 线 越 长 ， 


18 ”反之 则 越 短 。 三 组 试验 中 扩散 系数 从 Sm2.s 1 增加 10 倍 到 50 ms 4， 对 应 等 值 线 长 度 约 减 小 


19 六。 值得 注意 的 是 , 不 同 扩散 系数 并 没有 影响 相当 长 度 与 实际 长 度 间 的 比值 ， 两 者 比例 基 


20 ”本 保持 在 1.6 左右 。 


21 上 上述 分 析 表 明 ， 等 值 线 在 地 转 满 流 的 搅拌 作用 下， 充分 延展 伸 长 ,产生 多 尺度 精细 结构 
22 ”的 过 程 ， 同时 ， 小 尺度 扩散 不 断 地 抹 平 精细 结构 ， 抑 制 等 值 线 增长 ， 让 搅拌 作用 无 法 持续 产 


23 ” 生 精 细 结 构 。 两 者 在 清流 充分 发 展 后 达到 一 个 动态 平衡 。 在 这 个 平衡 下 ， 等 值 线 长 度 、 锋 面 


24 ”和 涡 丝 ( 即 梯度 因子 )、 混 合 效 率 等 特征 都 将 达到 一 个 稳定 值 。 这 一 点 从 等 值 线 的 时 间 演 变 


25 ”可 以 看 得 更 清楚 。 图 Sa 和 Sb 给 出 了 等 值 线 相当 长 度 和 实际 长 度 随 时 间 的 演变 情况 。 从 图 中 


26 ”可 以 看 到 , 两 个 长 度 都 有 约 100 天 左右 的 适应 和 调整 时 间 , 这 一 点 与 Shuckburgh 等 (2009) 
13 
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17 


的 结果 非常 一 致 。 这 段 时 间 内 两 个 特征 长 度 都 迅速 增加 ， 之 后 基本 达到 准 稳定 态 ， 演 变 趋 于 


缓慢 ， 说 明 大 约 100 天 后 锋面 和 涡 丝 在 小 尺度 扩散 的 作用 下 趋 于 饱和 ， 增 长 速度 受到 限制 。 


= 


但 是 ， 这 两 个 长 度 之 间 的 比值 “图 Sc) 则 基本 维持 稳定 ， 平 均值 大 约 在 1.6 到 1.7 左右 ,与 


前 面 的 分 析 一 致 。 


40 


20 


相当 纬度 (bea) 


-20 


-40 


-60 


0 100 200 300 0 100 200 300 0 100 200 300 
时 间 / 天 时 间 / 天 时 间 / 天 


图 5 (a) 标准 化 的 等 值 线 相当 长 度 ，'\b) 标准 化 的 等 值 线 实际 长 度 ，(c) 相当 长 度 与 实际 长 度 的 比值 。 
Fig. 5 Temporal evolution of (a) equivalent length Leg normalized by Lmin, (b) perimeter length L normalized by 


Lmin» and (c) ratio of Leg to L. 


2.2 等 值 线 长 度 与 混合 效率 、 扩 散 通 量 的 关系 


Haynes 和 Shuckburgh (2000) 给 出 了 等 值 线 三 个 特征 长 度 之 间 的 关系 式 : 


L, = (Val dj( 中 wela 


Ee (9) 
-7 (ival Val) = La: 


其 中 横 线 4=( 员 4dlj/( 种 加 为 沿 等 信 线 平均 算 子 ， 与 


(A), =(6,._, Al Val" dt) (G, [VaT a) 相似， 区别 在 于 是 否 用 两 根 等 信 线 间 的 宽度 


(Val) 加 权 。 对 上 式 除 以 及,, ， 可 以 得 到 混合 效率 的 表达 式 : 


E L aac 
Ma 区 vq)>1 (10) 
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从 (10) 式 不 难看 出 , 混合 效率 Mrz 包含 了 两 方面 的 贡献 , ESA SE PR EL 的 增长 ( 即 


性 / 尼 )， 二 是 示 踪 物 梯度 的 增强 〈 即 Val [Val )。 根 据 Cauchy-Schwarz 不 等 式 有 
Val ve|>1， 其 中 当 梯度 不 随 等 值 线 变化 时 《〈 即 沿 等 值 线 均匀 分 布 ， 如 图 la)， 等 号 成 立 


(Haynes and Shuckburgh, 2000)。 当 汕 流 搅拌 生成 锋面 和 涡 丝 时 ， 梯 度 沿 等 值 线 极 不 均匀 ， 


梯度 因子 |Vq| |Yd| 大 于 1， 这 就 是 相当 长 度 大 于 实际 长 度 的 原因 。 


相当 纬度 (ea) 


0 100 200 300 0 100 200 300 0 100 200 300 
时 间 / R 时 间 / 天 时 间 / 天 


图 6 (a) 混合 效率 Mp，(b) 等 值 线 增长 因子 L2/12in，(c) 梯度 因子 |vellval-1。 
Fig. 6 Temporal evolutions of (a) mixing efficiency Mg (b) perimeter length factor L2/L2,in, and (c) gradient factor 


|Vq||Vq|7?. 


图 6 给 出 了 混合 效率 及 其 两 个 贡献 因子 随时 间 的 演变 图 。 混 合 效率 , 也 就 是 有 效 扩散 与 


小 尺度 扩散 的 比值 Ken / K,， 代 表 了 在 涡 旋 搅拌 下 ,小 尺度 扩散 系数 被 放大 的 倍数 。 从 图 6a 


可 以 看 到 ， 在 地 转 滑 流 搅拌 下 ， 混 合 效率 可 以 提升 约 2-3 个 数量 级 。 三 个 混合 效率 高 值 区 分 


别 对 应 南北 半球 副热带 地 区 和 赤道 地 区 , 特别 在 北 亦 道 地 区 , 甚至 可 以 达到 1500 倍 。 不 过 ， 


赤道 流 场 经 过 了 散 度 订 ， 改 变 较 为 明显 ， 特 别 是 EKE 有 一 定 程 度 的 减弱 (图 略 )， 因 此 可 能 


还 低估 了 赤道 流 场 的 实际 混合 效率 。 另 一 方面 , 混合 效率 低 值 区 主要 在 南大 洋 地 区 (~50?S )， 


说 明 南 极 绕 极 流 是 混合 的 “障碍 ” 对 经 向 输 运 有 抑制 作用 (Ferrari and Nikurashin, 2010). 


对 比 图 6a 和 6b 不 难看 出 ， 混 合 效率 的 增加 ， 主 要 是 等 值 线 长 度 的 拉 伸 和 增长 ， 因 为 两 者 的 


15 


20 


21 


空间 分 布 几乎 完全 一 致 ， 大 小 也 在 2 个 数量 级 以 上 。 相 比 之 下 ， 梯 度 的 贡献 相对 较 小 (图 


6c)， 并 且 也 比较 固定 ， 大 约 在 2-3 倍 之 间 ， 而 且 spinup 的 时 间 更 短 ， 约 在 10 ALA, “DR 


度 锋 面 和 涡 丝 的 梯度 就 已 经 达到 饱和 并 维持 稳定 状态 。 


对 示 踪 物 控 制 方程 (1) 进行 纬 向 平均 ， 可 以 得 到 如 下 输 运 方程 : 


_ ] o _ 
oY = g cos g| -V Pe l (11) 
ôt cosy að aog j | 


方程 C1) 表明 示 踪 物 的 纬 向 分 布 受 等 号 右边 三 项 的 影响 ， 分 别 是 平均 流 输 运 、 涡 致 输 运 、 


以 及 小 尺度 扩散 带 来 的 输 运 。 通 常 进 一 步 化 简 的 方法 是 ,对 第 二 部 分 涡 致 输 运 利用 “ 顺 梯度 ” 


输 运 理论 进行 参数 化 得 到 : 


I E coso 5 F+ (Kug +r, ) 2 3 (12) 
ôt coso að 


RP K ay 是 经 验 性 的 涡 扩 散 系数 。 这 种 参数 化 方案 的 适用 性 可 以 通过 计算 有 效 扩散 理论 中 


的 扩散 通 量 来 验证 。 首 选 ， 本文 先 显 式 计算 了 涡 通 量 Y9 ， 然 后 基于 有 效 扩散 理论 计算 了 罕 


过 等 值 线 的 扩散 通 量 K gôg / (abpa) 。 这 个 扩散 通 量 与 参数 化 的 扩散 通 量 天 sw67 /(a09) 


9 着 相同 的 表达 式 , 但 两 者 物理 意义 不 一 样 。 前 者 是 严格 垂直 于 等 值 线 的 (类 似 于 跨 等 密度 


面 的 ) 通 量 ,， 且 有 效 扩散 KK 可 以 根据 等 值 线 的 拓扑 情况 准确 计算 出 来 , 因此 是 扩散 通 量 的 


准确 定义 (Nakamura, 1996); 后 者 是 严格 垂直 于 纬度 的 通 量 , 常常 因为 有 逆 梯 度 输 运 (如 Lu 


et al.,2016)， 导 致 天 saay 为 负 值 ， 即 使 在 顺 梯度 输 运 的 情况 下 ， 由 于 参数 化 的 高 度 不 确定 性 


而 使 Koy 无 法 准确 计算 。 正 因为 前 者 的 这 些 优点 ， 常 常 在 研究 中 用 来 栓 代 后 者 进行 分 析 


(Marshall et al., 2006; Thompson and Sallée, 2012 )。 


图 7 给 出 了 两 个 通 量 的 时 空 分 布 图 。 从 图 中 可 以 看 到 ， 涡 通 量 和 扩散 通 量 均 小 于 0， 表 


明 输 运 主要 从 北向 南 “ 顺 梯度 ”方向 进行 。 另 外 ， 涡 通 量 比 扩散 通 量 大 3 倍 左右 ， 并 且 两 者 


的 空间 分 布 也 有 明显 差异 首先 , 涡 通 量 大 值 中 心 在 亦 道 以 南 , 大 约 在 200 天 左右 达到 极 值 ， 
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而 扩散 通 量 中 心 在 (相当 ) 赤道 以 北 ， 且 一 直 在 随时 间 增加 。 虽 然 两 者 都 在 南北 半球 副热带 
各 有 一 条 大 值 带 ,但 它们 的 时 间 演 变 特征 相差 较 大 。 涡 通 量 演变 的 时 间 尺 度 较 短 , 瞬 变 性 强 ， 

这 也 使 得 9 (图 7a 等 值 线 ) 并 不 光滑 ， 相 比 之 下 ， 扩 散 通 量 (图 7b) 时 间 尺 度 较 长 ， 演 变 
缓慢 ， 并 且 基本 上 和 等 值 线 长 度 的 演变 〈 图 6b) 一 致 。 因 此 在 相当 纬度 坐标 中 (图 76)，g 
等 值 线 的 时 间 变 化 非常 平滑 。 更 重要 的 是 ， 涡 通 量 在 赤道 附近 和 北纬 40°N 附近 有 范围 不 大 
的 正 值 区 (图 7a 红色 等 值 线 )， 表 明了 局 地 的 逆 梯 度 输 运 〔 向 北 输 运 )。 因 此 ， 如 果 按照 经 
典 理论 将 涡 通 量 参 数 化 为 -V9 =K apg /(a6p) ,上 且 示 踪 物 的 南北 梯度 在 纬度 和 相当 纬度 
坐标 中 大 致 相同 〈 图 7 黑色 等 值 线 间隔 )， 不 难看 出 量 级 上 K yyy HOKU K p 的 3 倍 ， 并 


AK ay 还 会 出 现 局 部 的 负 值 ， 与 “ 顺 梯度 ”理论 相 了 矛盾。 


相当 纬度 (eq) 


o 50 100 


200 
时 间 / 天 时 间 / 天 


150 200 250 300 350 


7 基于 欧 拉 纬 向 平均 的 经 向 涡 动 通 量 和 沿 等 值 线 平均 的 扩散 通 量 对 比 图 
(a) BRAEG (105， 填 色 ) 和 58x (等 值 线 )， 其 中 [jx 表示 欧 拉 纬 向 平均 ，()' 表 示 对 应 的 距 平 ，(b) 基于 等 
值 线 平均 的 扩散 通 量 Kwag/agss/a 105, ME) 和 q(t Poa) 〔 等 值 线 )。 其 中 红色 等 值 线 为 通 量 的 0 等 值 线 。 


Fig. 7 Comparison of Eulerian zonal mean eddy flux and along-contour mean diffusive flux. 


(a) Eddy tracer flux v’q’* (10° color shadings) and a* (contours) based on Eulerian zonal mean (-)*; (b) diffusive flux 


Keff9q/9pea/a (10°) and q(t, Peq) based on the along-contour mean. Red lines indicate zero-flux contour. 


3 等 值 线 分 形 特征 
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(Kamenkovich et al., 2021; Qian etal.,2022 )， 第 二 类 则 从 谱 空 间 研 究 不 同 尺度 的 能 量 传递 和 


波 数 或 者 频率 的 变化 图 , 这 样 的 图 通常 使 用 对 数 图 , 图 中 曲线 的 斜 


串 级 规律 (Fox-Kemper and Menemenlis, 2008 )。 其 中 在 谱 空 间 研究 能 量 串 级 常常 给 出 能 量 随 


率 是 能 量 串 级 的 关键 参数 。 


有 趣 的 是 ,在 分 形 研究 中 ,分 形 维度 同样 也 是 定义 在 对 数 图 上 的 曲线 斜率 , 但 图 中 曲线 是 等 


值 线 长 度 随 不 同 网 格 尺 度 的 变化 。 两 者 的 相似 性 驱使 我 们 进一步 分 析 等 值 线 的 分 形 维度 , 并 


探讨 其 与 消 流 混合 的 关系 。 


本 文 以 0.1” 为 基准 ,采用 2" 倍 网 格 距 作为 “尺子 ”衡量 等 


值 线 的 长 度 ， 其 中 HX 0-7 


k8 组 结果 。 图 8 挑选 了 四 个 时 刻 等 值 线 的 分 布 图 ， 并 用 不 同 “ 尺 子 ” 作 为 度量 单位 计算 了 


等 值 线 的 长 度 。 从 图 中 可 以 看 到 ， 对 于 初始 时 刻 ， 不 同 “ 尺 子 ” 衡量 的 结果 基本 一 致 ， 说 明 


等 值 线 几 乎 没有 细节 信息 ， 因 此 等 值 线 是 简单 拓扑 曲线 ， 即 一 维 图 形 ; 随 着 时 间 的 增加 ， 等 


值 线 细节 信息 逐渐 丰富 ， 不 同 “ 尺 子 ” 衡 量 得 到 的 差异 也 就 逐渐 增 大 。 这 说 明 等 值 线 存在 尺 


度 不 同 的 长 度 信息 ， 逐 渐 显 现 分 形 特征 。 从 图 8d, 8f 8h 不 难看 出 ， 在 地 转 消 流 搅拌 下 的 


等 值 线 图 形 ,“ 斥 子 ” 越 大 ， 忽 略 的 细节 越 多 ， 得 到 的 长 度 估 计 值 就 越 得。 在 搅拌 作用 刚 开 


始 不 久 , 不 同 尺 度 间 计算 结果 的 差异 比较 均匀 ; 但 到 了 后 期 , 1 AF 


2 倍 网 格 距 的 差异 最 大 ， 


表明 大 量 细节 存在 于 2 倍 网 格 距 之 间 。 由 于 数据 原始 分 辨 率 的 限 


小 “尺子 ”衡量 的 结果 。 
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相当 纬度 (ea) 
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相当 纬度 (eq) 


0° 60°E 120° 180° 120°W 60° 示 踪 物 浓度 


相当 纬度 (eq) 


8 E) 示 踪 物 浓度 (无 量 纲 ) KESTER A) 等 值 线 长 度 〔 以 Lmin 为 单位 )。 审 图 号 为 
GS(2016)1665 

第 一 行 到 第 四 列 分 别 是 积分 1、20、100、365 天 后 的 结果 。 等 值 线 长 度 采用 了 “ 数 盒子 ” 方法， 展示 了 8 个 网 格 距 做 
尺子 度量 得 到 的 结果 。 


Fig. 8 (left) tracer distributions at different timesteps and (right) contour lengths normalized by Lmin- 


Rows one to four show the results of integration after 1, 20, 100, and 365 days. Contour lengths are calculated 


using ‘box-counting’ method with 8 grid boxes. 
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10|(a) 第 1 天 d [四 第 10 天 | Lo 第 50 天 


: 斜率 = 1.0519 â 斜率 = 1.4333 


2 1 1 1 1 L 


10 La) 第 100 天 


伸 率 = 1.4917 


log (等 值 线 长 度 / 网 格 距 ) 


1. 4887 ” 斜率 = 1.4829 


4 -3 -2 -1 0 -4 -3 -2 
-log (网 格 距 ) 
图 9 相当 纬度 53.3°S 对 应 的 等 值 线 长 度 Cm) 与 分 辨 率 Cm) 的 对 数 关系 

从 (a) 到 (i) 分 别 为 第 0、10、50、100、150、200、250、300、350 天 的 结果 ， 其 中 拟 合 斜 率 标注 在 各 个 子 图 内 ， 紫 色 阴 影 
为 95% 拟 合 置 信 区 间 ， 两 条 虚线 分 别 是 研究 时 段 内 最 大 和 最 小 拟 合 斜 率 的 曲线 。 


Fig. 9 Contour length (m) against resolution (m) on a log-log plot for the equivalent latitude of 53.3°S. 


Panels (a) to (i) are results at 0, 10, 50, 100, 150, 200, 250, 300, and 350 days. Linear-fit slopes are indicated in each panel. Purple 
shadings indicate 95% confidence intervals for the linear fit. Two dash lines are the fitted lines with maximum and minimum slopes 


during the time range of interest. 


有 了 不 同 分 辩 率 计算 的 结果 ， 可 以 用 log-log 图 和 线性 拟 合 来 估算 等 值 线 的 分 形 维度 。 


图 9 挑选 了 相当 纬度 53.3°S 对 应 的 等 值 线 作为 例子 来 展示 其 分 形 维度 的 计算 ， 因 为 在 这 个 


纬度 上 陆地 影响 最 少 〈 图 3b )。 由 于 选择 了 8 组 分 辨 率 ， 因 此 在 每 个 子 图 上 有 8 个 等 间距 样 


本 供 拟 合 趋势 。 从 图 9 不 难看 到 ,不 同时 刻 ， 等 值 线 长 度 与 分 辩 率 在 对 数 图 上 基本 满足 线性 


分 布 特征 。 特 别 地 ， 在 初始 时 刻 ， 由 于 等 值 线 为 直线 ， 因 此 斜率 、 即 分 形 维度 为 1 。 随 着 满 


流 搅拌 作用 的 产生 ， 和 斜率 逐渐 增加 ， 大 约 在 100 天 后 ， 达 到 一 个 相对 稳定 的 值 1.48， 并 有 小 


幅度 的 变化 。 虽 然 这 一 分 形 维度 离 二 维 图 形 还 有 较 大 差距 ， 但 是 比 经 典 的 Koch 雪花 曲线 的 


分 形 维度 1.26 要 大 ， 较 为 接近 Sierpinski 铁 热 的 分 形 维度 1.59。 
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Fig. 10 Spatio-temporal variations of the “box-counting” fractal dimension. 


以 上 是 用 一 个 相当 纬度 做 例子 。 实际 上 可 以 对 所 有 纬度 和 时 间 进 行 同样 的 计算 从 而 得 


到 分 形 维度 的 时 空 变化 情况 。 从 图 10 可 以 看 到 ， 分 形 维度 经 过 约 3 个 月 的 spinup 达到 准 稳 


定 状 态 ， 与 等 值 线 长 度 的 演变 一 致 。 前 面 所 挑选 的 53.3°S 则 处 在 维度 较 小 的 纬度 带 。 在 南 


北半球 副热带 地 区 和 赤道 地 区 ， 分 形 维度 可 以 达到 1.6 以 上 ， 特 别 是 赤道 地 区 ， 在 最 后 几 十 


天 甚至 接近 1.7, 较 Sierpinski 铁 垫 更 接近 二 维 图 形 。 这 也 说 明 亦 道 地 区 的 等 值 线 具 有 更 加 丰 


富 的 细节 信息 。 此 外 , 与 图 5 进行 比较 不 难看 出 ， 分 形 维度 与 等 值 线 长 度 的 演变 特征 非常 类 


似 。 因 此 ， 分 形 维度 可 能 与 混合 效率 有 特定 的 统计 关系 。 本 文 将 尝试 建立 分 形 维度 与 混合 效 


率 之 间 的 一 个 经 验 关 系 。 


通过 将 两 者 画 在 散 点 图 上 不 难 初步 估计 两 者 之 间 存 在 指数 关系 , 因此 不 妨 设 两 者 满足 如 


下 关系 式 : 


log(M,,)=aD,,, +0. (13) 


为 了 确定 a, b 这 两 个 经 验 系数 ， 定 义 残 差 R= 》 | aD,,,+b 一 log (Mj) |> RAB 


法 拟 合 (13) 式 两 边 ， 使 得 残 差 R 最 小 。 实 际 计算 中 发 现 ,不同 (相当 ) 纬度 的 拟 合 结果 有 
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一 定 差异 ， 因 此 本 文 将 全 球 划分 为 8 个 纬度 带 ， 间 隔 15°. 


1000 F(a) 纬度 带 [-60, -45] F(b) 纬度 带 [-45, -30] F(c) 纬度 带 [-30, -15] F(d) 纬度 带 [-15, 0] 时 间 / 天 


100 | R = 3070.05 L R= 195.92 | R=220.11 L R= 429.19 


AE S 
+ A log(ME) = 7.93Dmin-6.65 A 

Ta . r 
a ie T T | 


Af 1000 F(e) 纬度 带 [0, 15] 


log(Me) = 10.18Dmin-9.42 log(Me) = 10.23Dmin-9.50 log(Me) = 10.18Dmin-9.42 


T T T j 200 


E(g) 纬度 带 [30, 45] E(h) 纬度 带 [45, 60] 


E(f) 纬度 带 [15, 30] 


100 | R=382.90 L R = 309.17 L R = 1346.88 


hee s 
ys E 
log(Me) = 10.17Dmin-9.44 i 7 log(ME) = 9.27Dmin-8.05 ee log(Me) = 9.32Dmin-8.13 
1 1 1 1 E 1 1 1 E 1 1 1 A 1 1 1 0 
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图 11 分 形 维度 与 混合 效率 的 散 点 关系 
从 (a) 到 Ch) 分 别 为 不 同 纬 圈 范 围 的 结果 。 颜 色 表 示 散 点 的 时 刻 《〈 天 )。 图 中 
合 的 线性 方程 分 别 标 注 在 各 个 子 图 上 。 


Fig. 11 Scatter plot of fractal dimension Dmin and mixing efficiency Mg. 


phol 


线 为 最 小 二 乘 拟 合 的 直线 ， 拟 合 残 差 R 以 及 拟 


Panels (a) to (i) are results over different latitude bands. Colors indicates the time (day) relative to the first day. The black lines indicate 


a least-square linear fit, with residual R and the fitted equation labeled in each panel. 


图 11 给 出 了 8 SARE. TERRA BITE AFE, BU SE Reo 1 FE ALY 


左下 角 《〈 散 点 颜色 ); 经 过 一 年 左右 的 充分 搅拌 后 ， 散 点 基本 上 位 于 右上 角 。 从 散 点 分 布 趋 


势 和 拟 合 直 线 来 看 , 分 形 维度 和 混合 效率 的 对 数 基 本 满足 线性 关系 , 特别 是 在 消 流 充分 发 展 


后 。 但 是 在 发 展 初期 ， 特 别 是 分 形 维度 小 于 1.2 时 ， 散 点 和 拟 合 线 有 一 定 偏 离 。 此 外 ， 从 拟 


合 残 差 可 以 看 到 ， 在 45°S 到 30°N 范围 之 间 (图 11b-f， 残 差 都 小 于 500， 而 在 靠近 极地 的 


另外 三 个 纬度 带 〈 图 lla, g h), 拟 合 残 差 都 在 1000 以 上 。 特别 是 南半球 高 纬 地 区 (图 11a) 


的 残 差 可 达 3000 以 上 ， 不 少 散 点 显著 偏离 拟 合 线 ， 但 仍然 显示 出 线性 趋势 。 进 一 步 分 析 发 


现 这 些 条 带 状 的 离 群 值 〈outlier〉 随 纬度 增加 而 逐渐 远离 拟 合 线 ， 因 此 很 可 能 是 受到 地 形 和 


海 冰 的 影响 ， 拟 合 系数 a、4。 在 高 纬 地 区 有 显著 的 随 纬度 变化 的 特征 ( 表 1)。 不 过 ， 从 其 他 


纬度 带 相对 稳定 的 结果 (45”S~45”N) 可 以 得 到 a、5 在 最 小 二 乘 意义 下 的 估计 值 为 9.98 


22 


18 


19 


+0.006 和 -9.15 士 0.008， 即 : 


M , =exp(9.98D 


min 


-9.15), (14) 


并 且 该 关系 式 在 满 流 充分 发 展 起 来 后 吻合 度 更 高 。 这 一 点 在 线性 坐标 图 上 看 得 更 清楚 (图 
HE )。 


表 1 图 11 中 不 同 纬度 拟 合 系数 及 其 标准 误差 
Tab.1 Linear fitted coefficients and their standard errors in figure 11 for different latitude bands 


60°S 45°S 30°S 15°S 0° 15°N 30°N 45°N 

45°S 30°S 15°S 0° 15°N 30°N 45°N 60°N 

a 7.93 10.18 10.23 10.18 10.17 10.09 9.27 9.32 
+ an $ + are + 村 È + 

stderr 0.05 0.01 0.01 0.01 0.01 0.01 0.02 0.03 

b -6.65 -9.42 -9.50 -9.42 -9.44 -9.42 -8.05 -8.13 
土 ae al + 土 二 在 al + 

stderr 0.07 0.01 0.01 0.02 0.02 0.02 0.04 0.04 


R 3070.05 195.92 220.11 429.19 382.90 309.17 1346.88 1759.98 


HE: ay DAHA 3) 式 中 的 拟 合 斜 率 和 截 距 ，R 为 总 残 差 ，stderr 为 标准 误差 


Notes: a, b: fitted slope and intercept in Eq. (13), R: total residual of the fit, stderr is standard error 


本 文 利用 全 球 海 表 地 转 消 流 数据 驱动 示 踪 物 ， 考 察 示 踪 物 等 值 线 的 演变 特征 ,并 探讨 了 


其 与 混合 效率 、 扩 散 通 量 的 关系 。 研 究 发 现 ， 等 值 线 在 添 流 的 搅拌 下 会 产生 多 尺度 的 拉 伸 、 


扭曲 、 变 形 、 折 车 等 过 程 ， 这 些 过程 伴 随 着 丰富 且 精 细 的 涡 丝 和 锋面 结构 的 生成 。 在 最 初 的 


三 个 月 内 等 值 线 迅速 增长 , 之 后 由 于 小 尺度 扩散 的 平滑 作用 ， 等 值 线 的 增长 受 限 ， 并 在 搅拌 


和 扩散 两 个 动力 过 程 的 竞争 下 达到 动态 平衡 。 由 于 此 时 等 值 线 增长 10-20 倍 ,不同 性质 的 流 


体 接触 面 大 大 提高 ， 加 上 涡 丝 和 锋面 附近 梯度 的 显赫 增强 ， 混 合 效率 可 以 增 大 3 个 数量 级 ， 


即 涡 扩 散 系 数 可 以 比 小 尺度 扩散 系数 大 3 个 数量 级 。 定 量 计算 表明 , 等 值 线 的 增长 是 主要 的 


素 ， 梯 度 的 增强 则 是 次 要 因素 。 


| 
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此 外 , 本 研究 通过 计算 涡 动 经 向 通 量 和 有 效 扩散 通 量 ,发现 两 者 在 大 小 和 空间 分 布 型 上 


相差 较 远 ， 涡 动 经 向 通 量 甚 至 出 现 “ 上 梯度 ”输送 ， 因 此 传统 的 “ 顺 梯 度 ” 参 数 化 方案 并 不 


适合 。 


im 


得 一 提 的 是 ， 本 文 使 用 的 正 压 水 平 运动 经 过 纬 向 平均 后 , 涡 通 量 的 旋转 分 量 已 经 被 


H, 只 剩 下 辐 散 分 量 。 所 以 ,“ 顺 梯度 ”参数 化 方案 也 不 适合 涡 通 量 的 辐 散 分 量 , 与 Roberts 


和 Marshall (2000) 的 结论 不 一 致 。 因 此 ,“ 顺 梯度 ”参数 化 方案 的 合理 性 仍然 有 待 探 讨 。 


H FEER KE mE KEHE FER, 因此 等 值 线 具有 明显 的 分 形 特 


征 。 通 过 进一步 计算 等 值 线 的 分 形 维度 ， 发 现 其 时 空 演变 特征 与 等 值 线 长 度 、 混 合 效率 非常 


相似 。 因 此 ， 本 文 利用 最 小 二 乘 拟 合法 ， 得 到 了 分 形 维度 和 混合 效率 的 经 验 关 系 式 〈14 )。 


该 关系 式 在 河流 充分 发 展 后 吻合 程度 更 高 。 由 于 分 形 维度 可 以 从 遥感 观测 如 海 温 数据 等 ) 


的 等 值 线 直接 估算 得 到 ， 因 此 ， 该 经 验 关系 式 可 以 直接 推导 


混合 效率 ， 为 不 同 示 踪 物 的 混 


CE 


合 扩散 系数 提供 一 种 新 的 参数 化 思路 。 


由 于 本 研究 所 用 原始 数据 分 辨 率 在 0.25°*， 和 纬 圈 长 度 360° 之 间 只 跨越 了 3 个 数量 级 ， 


因此 其 范围 对 于 研究 能 量 串 级 和 等 值 线 分 形 维度 的 关系 仍然 不 够 宽广 , 无 法 考察 更 小 分 辩 率 


的 精细 特征 。 我 们 期 望 未 来 的 SWOT (surface water and ocean topography) 卫星 能 够 提供 更 


高 分 辨 率 的 海 表 流 场 数 据 ， 从 而 使 得 本 研究 能 够 延展 到 更 精细 的 尺度 。 在 更 高 分 辨 率 的 数据 


支持 下 (直至 分 子 尺 度 )， 不 难 想象 ， 小 尺度 扩散 系数 逐渐 接近 分 子 扩 散 ， 会 有 更 丰富 和 精 


细 的 涡 丝 和 锋面 ,等 值 线 在 搅拌 作用 下 会 更 加 长 ， 从 而 使 得 混合 效率 高 于 3 个 数量 级 。 但 希 


BARRA OH 在 更 精细 的 尺度 上 也 能 够 适用 。 


分 形 维度 是 数学 研究 的 前 沿 领域 , 尽 


mè 


有 研究 对 潮流 与 分 形 的 关系 进行 了 研究 , 但 两 者 


的 本 质 关 系 仍然 没有 理 


青 。 例 如 虽然 Navier-Stokes (N-S) 方程 与 分 形 有 关 ， 但 是 目前 还 没 


二 人 


TM N-S 方程 这 样 的 基本 的 运动 方程 中 直接 导出 分 维 的 理论 〈 黄 真理 , 2000)。 期 望 本 文 的 


研究 结果 能 够 为 这 方面 的 理论 研究 和 探索 提供 一 些 有 用 的 线索 和 参考 。 
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